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Abstract. Radio emission lias been implemented in 
the simulation software AIRES by M. A. DuVernois, 
B. Cai and D. Kleckner to perform a full Monte 
Carlo simulation of the radio emission of cosmic 
ray extensive air shower. Different sets of showers 
have been simulated. We present here the dependency 
of the radio electric field with different parameters 
(primary energy, radial distance, arrival direction, 
etc). We stress the effect of the polarization and the 
arrival direction, which should have a fundamental 
effect on the interpretation of experimental data. 
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I. Introduction 

The radio detection of cosmic rays has been first 
explored in the late 60's and early 70's H]. Thanks to 
progress in electronics, we currently assist to a renewal 
of experimental work with CODALEMA in France |2|, 
LOPES in Germany |3] as well as radio prototypes at 
the Pierre Auger Observatory, Argentina [4], |5|. 

Various models of emission have been proposed, 
and several are still under development such as the 
Coulombian emission of a charge excess, the transverse 
current or the geosynchrotron emission |I6|, Q, fS), ||9l , 
ifTOl . ifm . If a geomagnetic effect has been observed 
experimentally fT], f3|, |12|, a complete model which 
reproduces correctly all the observations is still to be 
developed. 

II. Radio emission AIRES 

The general form of the electric field produced by an 
accelerated relativistic particle is given by: 
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where v represents the velocity of the particle, n is the 
direction and R the distance from the emission point to 
the observer. The synchrotron emission corresponds to 
an acceleration induced by a magnetic field v — — v x 
B. 

Radio emission has been implemented in the sim- 
ulation software AIRES |13| (latter called ReAIRES) 
by M. A. DuVernois, B. Cai and D. Kleckner It 
is based on the Monte Carlo program AIRES in order 



to provide a realistic extensive air shower (EAS). For 
every electron and positron simulated, the electric field 
is computed based on the radiative term of Eq. [T] using 
the geomagnetic field at the location of the simulation. 
The refractive index of air n is approximated to unity. 
The static term of emission has also been implemented 
and has been checked to be negligible. 

The AIRES thinning algorithm is used to reduce the 
number of particles to simulate. The minimum number 
of particles to simulate mainly depends on the distance 
of the observation point to the shower core and on 
the considered frequency. The thinning level has been 
chosen appropriately for the following simulations. 

III. Shower SIMULATION 

For this analysis, showers were simulated on a virtual 
array of 40 antennas located up to 250 m around the 
shower core. The geomagnetic field is taken to Western 
Europe's values. The North-South (NS) axis is defined 
positive to the North, the East- West (EW) axis positive 
to the West and the vertical axis (VT) positive upwards. 
The shower direction is given by its direction of origin, 
with a zenith angle 9 and an azimuth angle (f) (0° = N 
and 90° = W). During the analysis, the radio signals 
are filtered in the 23-83 MHz band, even though the 
following results do not depend much on this band. We 
latter call signal the amplitude of the peak of the filtered 
signal, which can be positive or negative. 

A typical radio footprint (10^^ eV proton induced 
EAS coming from 30° zenith and 45° azimuth, i.e. 
North- West) is given on Fig. [T] 

The electric field is generally axially symmetric 
around the shower axis. This is visible on Fig. |2l where 
the signals from the same shower are plotted against the 
radial distance to the shower axis in the shower plane. 
The signals are fitted with an exponential function of the 
form: 

E^Ea cxp (2) 

where Eq represents the field on the axis, R is the 
radial distance and i?o is the characteristic decreasing 
distance of the field. Eq and Rq are adjusted during the 
fit. For most of the simulated showers, the signals are 
well describes by the exponential function. The i?o in 
the three polarizations are very similar. 

For particular arrival directions of the shower, the 
footprint in a given polarization can exhibit an asym- 
metric pattern. A vertical shower for instance presents 
two changes of sign for the NS field, and one change of 
sign for the vertical field. However the EW field, which 
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Fig. 1. Footprint of the electric field, in each polarization, for a lO^'^ eV shower coming from 9 = 30° and = 45°. The dots indicate the 
points where the field is actually calculated. We can note a change of sign between the polarizations. 
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Fig. 2. Lateral distributions of the signals of Fig. [T] fitted with exponential functions. 



value is typically ten times greater than the other fields 
here, is axially symmetric; as illustrated on Fig. [3] More 
generally, these kinds of patterns only appear when the 
signal in this polarization is very weak compared to the 
other polarizations, and can thus be ignored in a first 
order description of the electric field. It actually happens 
only when the projection of v x B in the considered 
polarization is very small, as discussed latter 

IV. Arrival direction and v x B polarization 

In order to study the effect of the arrival direction on 
the radio signals, sets of 121 showers of fixed energy 
and coming from various directions (from 9 = to 
75°) have been simulated. The amplitudes Eq in all 
three polarizations have been extracted for each shower 
Their values present strong variations upon the arrival 
direction. This is shown for 10^^ eV proton primaries as 
a function of the arrival direction on Fig. 2] The NS and 
VT terms are very similafl there are null for showers 
coming from North and South, are positive for showers 
coming from East and negative for showers coming from 
West. For the EW term, it is positive for showers coming 
from a larger North area, null on an East-West curve 
containing the geomagnetic field direction and negative 
on South. 

'There is actually just a tangent of the geomagnetic inclination 
between the two. 



This is simply explained by a linear dependence 
between the electric field and the vector cross product 
of the direction of the shower v and the geomagnetic 
field direction B: E oc — v x B, as predicted by some 
analytical models. It is the case for instance for the 
synchrotron radiation when the electric field is observed 
near the particle motion axis. Thus the skymaps pre- 
sented in Fig.lHcan be parameterized simply considering 
this vector cross product and a single function of the 
zenith angle: 

E = -F(6l).vxB (3) 

The vector cross products represents the emission mech- 
anism, and the function F{0) contains the shower de- 
velopment, the longitudinal distance to the shower, the 
atmospheric depth, etc. This function F{6) has been 
adjusted using a function of the form A.{1 + B.6)/{1 + 
exp{^-^)). The functions obtained are shown on Fig.|6] 
for proton primaries at different energies ranging from 
10^^ to 10^^ eV. The F{9) functions are almost propor- 
tional to the primary energy. The remaining differences 
may be interpreted in term of shower development: 

> A low energy shower reaches its full development 
higher in the atmosphere, thus produces less elec- 
tric field on the ground. When the zenith angle 
increases, the shower gets even farther, and the field 
is weakened. 

> A vertical high energy shower hits the ground 
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Fig. 3. Same as Fig. \l\ but for a vertical shower. 




Fig. 4. Skymaps of Eg from ReAIRES: amplitude Eo{6, 4>) as a function of the anival direction of the simulated showers (indicated with 
dots), in the NS (left), EW (center) and VT (right) polarization. 
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Fig. 5. Skymaps of F{9).{v X B), where F{9) is adjusted to fit the skymap of Fig.|4] 



before it reaches its full development. Increasing 
the zenith angle will allow the shower to develop 
above the ground and increase the electric field on 
the ground. After a given zenith angle, the shower 
will become too far from the observation point and 
the signal will start to decrease again, as for a low 
energy shower. 



This corresponds to the tendency observed in Fig. |6l the 
global maximum tends to be reached at higher zenith 
angle as the energy increases. 



V. Overall parameterization 

A. Energy dependency 

As already briefly discussed, the signal amplitude is 
quasi proportional to the energy of the primary cosmic 
ray. The simulation of a set of 41 showers in a fixed 
direction {0 — 30°, = 45°) with energies ranging from 
10^^ to lO'^" eV gives the relation E cx Energy^-^^, thus 
we will simply consider the proportionality. 

B. Characteristic distance Rq{9) 

As the zenith angle 9 increases, the shower gets farther 
from the observation point and the signal lateral distri- 
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Fig. 6. F{6) for proton primaries at different energies. Tlie reference 
function is tlie red full line, at lO'^'^ eV. The other functions have 
been rescaled by 10^ at 10^^ eV (mixed black curve), by 10~^ at 
10^® eV (dashed blue curve) and by 10~* at 10^^ eV (long dashed 
green curve). 



bution on the ground is flattened, i.e. the characteristic 
distance Rq of Eq. |2] increases. This effect is added to 
the broadening of the footprint on the ground due to 
the projection effect. The distribution of Rq{9) for the 
10^^ eV showers used to produce the Fig. |4]is given in 
Fig. |2] The distribution is fitted here with the function 
Ro{e) = 61.7 (1 + 3.43 10-^6i3-22) m. 




Fig. 7. Ro{e) the events of Fig.g] 

These distances are sensibly smafler than the experi- 
mental measurements, as for most current models, and 
this point needs to be clarified. However, the dependence 
over the zenith angle still is interesting. 

C. Global formula 

Different dependencies of the electric field have al- 
ready been exposed: 

• proportionality to the primary cosmic ray energy 
due to the number of particles 



• proportionality to the — v x B linked to the emission 
mechanism 

> a dependence on the zenith angle via F{6), to 
interpret in terms of atmosphere and shower de- 
velopment 

• an exponential decrease with the distance of to the 
shower axis via Rq{9) 

To generalize we can add a proportionality to the am- 
plitude of the magnetic fielq^. This gives the following 
formula for the electric field produced by an HAS around 
10^'' eV in a bandwidth of Ai^ MHz: 

Ai^ \B\ Energy 
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VI. Discussion 



Synchrotron radio emission was implemented in 
AIRES to perform a full Monte Carlo simulation of the 
radio emission of cosmic rays EAS. We have detailed 
here an analysis and parameterization of the electric field 
simulated for different configurations. If some difference 
with experimental measurements still need to be under- 
stood (electric field too high, Rq too small), important 
characteristics of the electric field can be extracted (as 
the —V X B polarization) and can guide experimental 
analysis, as for instance in lfT4l . 
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